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Abstract 

We solve numerically the equations of motion for the collapse of a 
shell of baryonic matter falling into the central regions of a cluster of 
galaxies, taking into account of the presence of the substructure induc- 
ing dynamical friction. The evolution of the expansion parameter a(t) 
of the perturbation is calculated in spherical systems. The effect of 
dynamical friction is to reduce the binding radius and the toatal mass 
accreted by the central regions. Using a peak density profile given by 
Bardeen et al. (1986) we show how the binding radius of the perturba- 
tion is modified by dinamical friction. We show how dynamical friction 
modifies the collapse parameter of the perturbation slowing down the 
collapse. 
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1. Introduction 



In the most promising cosmological scenarios structure formation in 
universe is generated through the growth and collapse of primeval den- 
sity perturbations originated from quantum fluctuations (Guth & Pi 
1982; Hawking 1982; Starobinsky 1982; Bardeen at al. 1986) in an in- 
flationary phase of early universe, perturbations is due to gravitational 
instability. The statistics of density fluctuations originated in the in- 
flationary era arc Gaussian and it can be expressed entirely in terms 
of the power spectrum of density fluctuations: 

P{k) =< \S^\' > (1) 

where 

5k = J d^kexp{-ik:>c)5{:>c) (2) 

S{^) = f^^^ (3) 
Pb 

and Ph is the mean background density. In biased structure forma- 
tion theory it is assumed that cosmic structures of linear scale it! 
form around the peak of density field, 5(x), smoothed on the same 
scale. Density perturbations evolve towards non-linear regime bacause 
of gravitational instability and it breaks away from general expansion 
at: 

1/2 

(1 + zf^ (4) 
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where z is the rcdshift, 5 is the overdensity within r. When 5 ~ 1 
the density perturbation begins to recollapse. The collapse time, Tco, 
depends on the characteristic of initial fluctuation field as average over- 
density and on the environment in which the perturbation is embedded. 
This last features depends on the cosmological scenario . The mean 
characteristics of the structure of the structures which form around 
the peak of the density field depend on the spectrum, P{k), which 
in turn depends on the matter that dominates the universe (CDM, 
WDM, HDM). In this context the most succesful model is the biased 
CDM model (Liddlc & Lith 1994) based on a scale invariant spectrum 
of density fluctuations growing under gravitational instability and on 
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the hypothesis that the universe is dominated by cold dark matter. 
A very simple model for accretion of matter in cluster of galaxies was 
first investigated by Gunn & Gott (1972). In their accretion model the 
density profile is given by: 

Rf 

Pi = Pei + {Pci + P+- Pei)^ r > Ri 

pci + p+ n < Ri (5) 

where p^ and p^ are the critical and external density, p+ gives the 
density surplus present in the radius Ri respect to the critical density. 
The perturbation with p > p^ reach a maximum radius r^^ = | at a 
time: _ 

7^2 = ^(1±^ (6) 

' Hi (53/2 ^ ^ 

where 5 is the overdensity in the radius r and is the Hubble param- 
eter at time tj. After Tco/2 the matter collapse and its infall is radial. 
The perturbation needs a time T^o for the total collapse. Gunn & Gott 
model is a oversimplification of the perturbation evolution. In Gunn & 
Gott model there are no tidal interactions among the shell of baryonic 
matter and the external density perturbations and it is supposed that 
there is no substructure, (collapsed objects of lenght less than that of 
main perturbation). One of the features of CDM models is an abun- 
dant production of substructure. In this scenario structure formation 
goes from bottom to up through gravitational clustering, merging and 
violent relaxation of small scale substructure (White & Rees 1978). So 
the matter inside a given region is clumped in a hierarchy of objects 
of various dimensions. The time a given lenght goes nonlinear can be 
obtained from the condition on mass spectrum, (7o(M): 

ao{M) ^1^1 + Zni (7) 

where: 

ao{M) = ^J P{k)k^W{kRf)dk (8) 
beeing Rf the filtering scale and W{kRf) the window function: 
W{kRf) = 3[5m(A;i?/) - kRfCosjkRf)] 
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So a subgalactic scale collapse at a rcdshift: 

30 

1 + Znl^— (10) 

(Silk & Stebbins 1993), where b is the biasing parameter. In particular 
a perturbation of IO^^Mq collapses at a redshift z ~ 0.02 while pertur- 
bations in the range lO^M© — lO^M© collapse almost at the same 2; ~ 18 
because the mass variance in this region varies only of a factor 3 (Rees 
1986). A part of this last pcrurbations has a cross section too little for 
gravitational merging to be dcstroied (Rees 1986). Therefore a mass 
perturbation of lO^M© — lO^M© can survive untill the cluster enters 
nonlinear phase at z ~ 0.02. Because of the presence of substructure 
in a cluster Gunn & Gott model needs a revision. In fact substructure 
acts as a source of stochastic fluctuations in gravitational field and in- 
duces dynamical friction (Antonuccio & Colafrancesco 1994) and this 
produces a modification of the motion of shells of baryonic matter in 
a density perturbation. In particular dynamical friction delays the 
collapse of low density perturbations {S ~ 0.01) (Antonuccio & Co- 
lafrancesco 1994). The plan of the paper is as follows. In section 2 
by numerical integration of equations of motion of a shell of baryonic 
matter made of galaxies and substructrure of lO^M© — lO^M© we show 
how the expansion parameter of the same varies with time. In section 3 
we show how dynamical friction affects the binding radius of a cluster. 



2. Modification of the expansion parameter of a shell. 

The equation of motion of a shell of baryonic matter around a maximum 
of the density field, neglecting tidal interactions and substructure, can 
be expressed in the form: 

_ GM 

d^~~^ ^^^^ 

(Peebles 1980, eq. 19.9), where M is the mass enclosed in the proper 
radius r{t). Using Gunn & Gott notation the proper radius can be 
written as: 

r{ri,t) ^ a{ri,t)ri (12) 
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where is the initial radius, a(rj, t) is the expansion parameter of the 
sheU. At the initial time ti the initial condition is given by a(rj, tj) = 1. 
In presence of substructure it is necessary to modify the equation of 
motion Eq. (|TT]) because of the graininess of mass distribution in 
the system, due to substructure. In a material system gravitational 
field can be decomposed into an average field, Fo(r), generated from 
the smoothed out distribution of mass, and a stochastic component, 
Fstoch{r), generated from the fluctuations in number of the neighbour- 
ing particles. The stochastic component of the gravitational field is 
specified assigning a probability density, W{F), (Chandrasekhar & von 
Neumann 1942). In a infinite Homogeneous unclustered system iy(F) 
is given by Holtsmark distribution (Chandrasekhar & von Neumann 
1943) while in inhomogeneous and clustered systems iy(F) is given 
by Kandrup (1980) and Antonuccio & Barandela (1992) respectively. 
The stochastic force, Fgtoch, in a self-gravitating system modifyes the 
motion of particles as it is done by a frictional force. In fact a particle 
moving faster than its neighbours produces a deflection of their orbits 
in such a way that average density is greater in the direction opposite 
to that of motion causing a slowing down in its motion. The modified 
equation of motion of a shell of baryonic matter can be written in the 
form: 

d'r GM 

(Langevin 1902, Saslaw 1985) where t] is the coefficient of dynamical 
friction. Supposing that there are no correlations among random force 
and their derivatives we have: 

^ Jd^FWiF)F'TiF) 
' 2 < f 2 > ^ ^ 

(Kandrup 1980), where T{F) is the average duration of a random force 
impulse, W{F), is the probability distribution of stochastic force that 
for a clustered system is given by Antonuccio & Barandela (1992). 
Equation (|T3p in terms of a(rj,t) and 
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can be written as: 

(fa _ ATiGpcijl + 6i) da 

dt^ ~ a\t) dt ^ ' 

where pd is the background density a time tj and 5i is the overdensity 
within Tj. Using the parameter r = ^ the Eq. ([T6|) may be written 
in the form: _ 

d^a _ A-nGpcijl + ^i) ^^ da 

JT^ - ^) ~ "^^"'d^ ^^^^ 

Referring to the calculation of r], given by Antonuccio & Colafrancesco 
(1994) and using equation ( p!7D the time evolution of the expansion pa- 
rameter, a(r), of the shell can be obtained. Equation ( p!7[ ) can be solved 
looking for an asymptotic expansion. This was made by Antonuccio & 
Colafranceasco (1994). In their paper they gave only an expression for 
the collapse time Tc. Equation (|l3) can be also solved numerically. To 
this aim we used a Runge-Kutta integrator of 4*^^ order. We studied 
the motion of a shell of low density, 5 = 0.01, typical of a perturbatoon 
present in the outskirts of a cluster of galaxies. The initial conditions 
were chosen remembering that in expansion phase the shell moves with 
Hubble flow and that at the maximum of expansion the initial velocity 
is zero. In figure 1, we show the expansion parameter a(r) versus r 
both when the dynamical friction is taken into account and when it is 
absent. As we can see dynamical friction slows down the collapse of 
the shell of matter in agreement with the analytic calculation of the 
collapse time, Tc, of a shell in which dynamical friction effect is taken 
into account given by Antonuccio & Colafrancesco (1994). 

3. Binding radius of a cluster in presence of dynamical 
friction. 

In biased galaxy formation theory structures form around the local 
maxima of the density field. Every density peak binds a mass m that 
can be calculated when we know the binding radius of the density peak. 
The radius of the bound region for a chosen density profile 6{r) may 
be obtained in several ways. A first criterion is statistic. The binding 
radius of the region, r^, is given by the solution of the equation: 

< S{r) >=< {6- < 6 >)2 >i/2 (18) 
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(Ryden 1988). At radius r « the motion of particles is predeminant 
toward the peak while when r » the particle is not bound to the 
peak. Another criterion that can be used is dynamical. It supposes 
that the binding radius is given by the condition that a shell collapse 
in a time, smaller than the age of the universe to- 

Te(r) < to (19) 

(Hoffmann & Shaham 1985). This last criterion, differently from 
the previous one, contains some prescriptions particularly connected 
with the physics of the collapse process of a shell. For this reason we 
used it to calculate the binding radius. The time of collapse, Tc(r), at 
radius r can be obtained solving numerically equation ( [TE| ) for different 
values of 5, from a given density profile 5{r). We use the average density 
profile given by Bardeen et al. (1986): 

z/^(r) ^(z^7,7) 



6{r) = A 



e(0)i/2 ^^(0)1/2(1-^: 



2^ 



(20) 



where A is a constant given by the normalization of the perturbation 
spectrum, P{k), v = ^(r) is the correlation function of two 

points, 7 and i?* two constants obtainable from the spectrum (see 
Bardeen et al. 1986) and finally 6'(7z/, 7) is a function given in the 
quoted paper (eq. 6.14). Given the average density profile the average 
density inside the radius r in a spherical perturbation is given by: 

5 = A r dx6{x)x^ (21) 



Jo 

We calculated the time of collapse, Tco(r), using Eq. (^ with the density 
profile given in Eq. (^). We repeated the calculation of Tco(r) for 
1.5 < u < 3 and we applied the condition given in Eq. ( |19D to the 
curves Tco(r) previously obtained. The result is the plot in Fig. 2 for 
the binding radius rf, versus z/. 

4. Conclusion 

In the first part of this paper we show how the dynamical friction 
slow down the collapse of a spherical density perturbation. The effect 
grows with increasing of the coefficient of dynamical friction rj. In the 
second part, we show how dynamical friction reduces the binding radius 
(see fig. 2). 
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Figure 1: Temporal evolution of the expansion parameter a(r) of a shell of 
matter made of galaxies and substructure versus r. The dashed line shows 
a(r) when dynamical friction is absent, while the solid line is the same when 
dynamical friction is taken into account. We assume a cluster radius of 
Rci = 5h~^Mpc, a central overdensity 5 = 0.01 and a total number of peaks 
of substructure Ntot = 10^. 
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Figure 2: Variation of the binding radius Vb with i/. The sohd Hne is the 
binding radius in absence of dynamical friction, while the dashed line is the 
same as in presence of dynamical friction. The filtering radius used is , 

Rf = 1 Mpc 



10 



